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Using recently developed eSicient symbolic manipulations tools, we present a general gauge- 
invariant formalism to study arbitrary radiative (/ > 2) second-order perturbations of a 
Schwarzschild black hole. In particular, we construct the second order Zerilli and Regge- Wheeler 
equations under the presence of any two first-order modes, reconstruct the perturbed metric in terms 
of the master scalars, and compute the radiated energy at null infinity. 

The results of this paper enable systematic studies of generic second order perturbations of the 
Schwarzschild spacetime. In particular, studies of mode-mode coupling and non-linear effects in 
gravitational radiation, the second-order stability of the Schwarzschild spacetime, or the geometry 
of the black hole horizon. 
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I. INTRODUCTION 

Nonlinearities are a characteristic feature of General 
Relativity. Classifying the possible types of nonlinear- 
ity can help us much in understanding and interpreting 
the results of observations or simulations of gravitational 
phenomena. For instance, the space-based gravitational- 
wave observatory LISA will be able to detect quasinor- 
mal mode coupling and frequency mixing in the ring- 
down phase of the collision of two supermassive black 
holes at distances of 1 Gpc. Perturbation theory offers a 
systematic study of nonlinearities, at least for moderate 
couplings, and it is hence worth developing a complete 
formalism for perturbations in General Relativity. 

During the last three years we have prepared a com- 
bination of ideas and tools on which such systematic ap- 
proach to metric perturbation theory in General Rela- 
tivity can be based. This includes, among other things, 
general explicit formulas for the perturbations of curva- 
ture tensors at any order [l| , a general analysis of gauge- 
invariance , and the construction of specialized and ef- 
ficient tensor computer algebra tools to handle the enor- 
mous equations of perturbative General Relativity [3, Q] ■ 
For definiteness, and because this is the most standard 
scenario in astrophysical applications, we have generally 
focused on spherical background spacetimes. 

The purpose of this article is to apply that general 
formalism to a Schwarzschild background. Linear per- 
turbations of Schwarzschild have been studied since the 
pioneering work of Regge and Wheeler Q , and later Zer- 



illi and Moncrief Q, in which they were already able 
to identify and describe the two polarizations of gravi- 
tational waves, giving two decoupled wave equations for 
them. 

The first studies of second-order black hole perturba- 
tions were carried out in the seventies [1, Q in order to 
study the non-linear stability of the Schwarzschild solu- 
tion. In the mid- nineties, motivated by the close- limit 
approximation gauge-invariant second-order gener- 
alizations of the Regge- Wheeler and Zerilli-Moncrief for- 
malisms were developed and applied to a variety of close 
limit-type initial data for binary black holes; see [ll|, [13 
and references therein. More recently, the relevance of 
second-order perturbation theory on an extreme mass- 
ratio inspiral (EMRI) has been evaluated [l^ . In general, 
those references have focused on the study of the self- 
coupling of a particular first-order mode and how that 
generates some second-order mode. 

We present here a complete, gauge-invariant Regge- 
Wheeler-Zerilli-Moncrief like formalism for arbitrary I > 
2 first and second-order perturbations of a Schwarzschild 
black hole. In particular, we derive the general first and 
second order Regge- Wheeler and Zerilli master variables 
and the equations they obey. We also reconstruct the 
perturbed metric in terms of those scalars, as well as 
computing the radiated energy at null infinity. 

Most of our results are not only gauge invariant but 
also covariant; that is, independent of the background co- 
ordinates. This is analogous to the covariant formalisms 
for linear perturbation theory in references [l5| and 
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0, [i3|- This is a necessary step to, for example, test 
the geometry near the horizon or null infinity, for which 
Schwarzschild coordinates are inadequate. Other results 
are still given in Schwarzschild coordinates, for simplicity. 

The organization of this paper is as follows: we start 
in Section |TT] by reviewing the formalism of [U, 0| • Sec- 
tion [nT] presents the high order Regge- Wheeler and Zer- 
illi equations and Section IIVI sketches how we compute 
the radiated energy. Setting up the arena for the sec- 
ond order treatment. Section IVl rederives in a compact 
yet complete way a gauge-invariant version of the Regge- 
Wheeler-Zerilli first order formalism. Finally, Section 
|Vl] presents our general second order Regge- Wheeler and 
Zcrilli equations, their sources, and the computation of 
the radiated energy in terms of our second order Regge- 
Wheeler and Zerilli functions. 



where gAB is a metric field and r a scalar field, both 
on the manifold Ai^, and 7ab is the unit metric on the 
sphere S*^ . The respective coordinate systems on the two- 
dimensional manifolds can still be freely specified. 

We define the following notation for the covariant 
derivatives associated with each metric: 



gAB\C = 0, lab-.c = . (2) 

For future reference, we define the vector va = r\Al'''- 



II. HIGH-ORDER GERLACH AND SENGUPTA 
FORMALISM 

A. Background spherical spacetime 

This section briefly summarizes the formalism intro- 
duced in references [ll, Q to deal with high-order per- 
turbations of a spherical spacetime. This formalism can 
be regarded as the generalization to higher orders of the 
Gerlach-Sengupta linear formalism [T3 Il5j|. 

We start by decomposing the spacetime manifold 
as a product M = x S*^, where is a two- 
dimensional Lorentzian manifold and S'^ the two-sphere. 
We will use Greek letters {n^v, . . .) for four-dimensional 
indices, capital Latin letters [A, B, . . .) for indices on Ai^ 
and lowercase Latin letters (a, 6, ... ) for indices on the 
sphere. With this notation, and choosing an adapted co- 
ordinated system, a general spherically symmetric back- 
ground metric is given by 

g^,^ix^,x'^)dx''dx'' = gAB{x'^)dx^dx^ (1) 
+ r^ix'')^abix'')dx''dx', 



B. Non-linear perturbations 

In perturbation theory one works with a family of 
spacetimes {M{s), g{e)) depending on a dimensionless 
parameter e. The background spacetime is the member 
of this family for which e = 0, and it is assumed to be a 
known solution of the Einstein equations. Performing a 
Taylor expansion around the background metric g^j^, 



n=l 

defines the perturbations ^"^h^^, which are tensors on the 
background manifold. In practice, we will always work 
at a given finite maximum order N , truncating the series 
at that order. 

Using the basis of tensor harmonics described in ap- 
pendix 13 we decompose the perturbations of the metric 
in the following way. 
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Perturbations of all relevant curvature tensors in term 
of metric perturbations have been given in [l[ , in partic- 
ular those of the Einstein tensor. The latter have also 
been decomposed there in spherical harmonics at second 
order, giving full equations of motion for the harmonic 
coefficients of the metric. 

Now we proceed to extract the gauge freedom present 
in the previous decomposition Per each perturbative 



order n, four of those ten components are free. 

The most natural gauge choice in spherical symmetry 
is the one introduced by Regge and Wheeler (RW) Q for 
linear perturbations. Here we use those same conditions 
at all orders, for I > 2: 

i-inrA = 0, <">Gr = 0, ^">/ir = 0. (5) 

Alternatively, it is possible to construct explicit gauge- 
invariant combinations of those harmonic coefficients, as 
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explained in 01- The idea is simply to select a gauge 
and define the gauge invariants to be the expression of 
the general gauge transformation of the metric functions 
from an arbitrary gauge to that selected gauge. Around 
spherical symmetry, the requirement that the gauge in- 
variants are local (almost) uniquely selects the RW gauge 
as preferred reference gauge, and this is what we will use. 

Schematically (see [3 for more details) the gauge in- 
variants associated the RW gauge are given by 
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(6) 

(7) 

(8) 



where the terms ^"^Rab, ^"'i?^ a-nd arc polynomial 
combinations of the lower order harmonic coefficients 



I A, 



'^GT, V}, with k < 



(9) 



all vanishing in RW gauge, by construction. Because 
of the k < n condition, those terms arc not present 
in linearized perturbations (n = 1). For second order 
perturbations (n = 2), they are quadratic in the first- 
order quantities {^'^HJ^a, ^'^G[", ^'^hY"}, and so forth 
for higher n. 

Once more, the values of ^"^ICab, *"^/C, in any 

gauge coincide with their values in the RW gauge. Be- 
cause of this we can do calculations in the RW gauge, 
yet still recover the form of any expression in a generic 
gauge by making use of the definitions ([SJH]). Wc will 
take advantage of this fact later. 



III. HIGH ORDER REGGE- WHEELER AND 
ZERILLI EQUATIONS 



equations to two wave equations for two TTiaster functions, 
one of odd / axial parity and another one even / polar par- 
ity, and this is true at all perturbative orders. These two 
functions fully describe the gravitational content of the 
system, and actually the whole metric perturbation at 
order n can be reconstructed from the master variables 
at orders k < n. 



A. The even/polar parity sector 

We define the nth order Zerilli function as the following 
combination of polar harmonic coefficients. 



3M + A 



(«^^">/CAi3- ^"%m)«^ + ^'"'/C, (11) 



where A = i(^— 1)(/4-2). This variable is given in terms of 
the nth order gauge invariants tied to the Regge- Wheeler 
gauge. To recover the form of the Zerilli function in a 
different gauge, we just need to replace the invariants 
by their explicit form in terms of a generic gauge ([6]- 
H]). Then, the Zerilli variable takes the following form in 
Schwarzschild coordinates: 



(2M - r) 
3M + Xr 



+ i;(/ + i)rf">G + 



(12) 



where the subindiccs t and r stand for components of the 
different tensors in these coordinates. In the previous ex- 
pression collects the terms { ^"^Ra, ^"^R}, 
and hence it is itself a polynomial in the lower order vari- 
ables Notice that in Eq. the last three terms, 
including are zero when imposing the RW gauge. 

The Zerilli master function satisfies the following wave 
equation 



{"}< 



(13) 



From now on we shall assume that our background so- 
lution is the Schwarzschild spacetime. The perturbative 
formalism summarized in the previous section allows us 
to work with arbitrary coordinates on the reduced man- 
ifold Ai^. However we will sometimes use Schwarzschild 
coordinates {t, r) as intermediate tools for computing ex- 
pressions which in the end will be valid in any asymptot- 
ically flat background coordinates. Then we will use the 
following shorthands for coordinate derivatives acting on 
any object uj: 



du! 
'dt' 



dr 



(10) 



When studying perturbations of the Schwarzschild 
spacetime, it is possible to reduce the perturbed Einstein 



The source term '"^5$ depends on the lower order per- 
turbations, while the potential is 



+ 6M r^X{X + 2) + 3Mir - M) 

= ~~;r2 ~?, (rA + 3Af)2 ■ 



(14) 



In particular, when using the tortoise coordinates 
(i, r*) [with r* ^ r + 2A/ln (r/ (2M) - 1)] the differential 
operator takes the following simple form 



2M 
r 



Q^2 Q^*2 



(15) 



Obviously, the first-order source ^^'S** is zero. The 
second order one can be given as 
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(16) 



where (Z, m) and {I, ffi) are a pair of two first-order modes 
which contribute to the second-order (/, m) mode. The 
sources S in the previous expression are expficitly given 
in reference Q. The polarity sign e is defined as e = 

Note finaUy that the definition of the high-order Zer- 
iUi function (fTTj) is essentially determined up to addition 
of low-order gauge-invariant terms. That is, the addi- 
tion of such low-order terms would keep the same form 
of the Zerilli equation ([13]), in particular with the same 
potential Vz, but would change the source '"'Sij,. The 
definition given in ((TT|) is just the simplest possibility, 
and follows We will later make use of this freedom. 



B. Odd/axial parity sector 

The Gerlach and Sengupta (GS) master scalar is de- 
fined as the rotational of the axial invariant vector Kji/r'^, 



{"} 



(17) 
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Like the Zerilli function, it is given in terms of some of the 
RW gauge-invariants ([6][8]). In a generic gauge it takes 
the form 



AB 



On, 



(18) 



\B 



where '"'Qn is a source term that depends on lower-order 
perturbations, and which is zero in the RW gauge. 
It obeys the GS master equation. 



_L(,..-.n)i.^ 



axial sector. As we will see, all metric components can 
be algebraically reconstructed from this scalar. 

We will use the above definition for ^"'11 for historical 
reasons. But, in fact, a better choice [2l| is the rescaled 
^"*n = *"*n, because its evolution equation has no 
first-order derivatives. 



(21) 



This equation is valid in any spherically symmetric back- 
ground and the potential is given by 

^RW = ^^-^(l-5^VnB)- (22) 



For the Schwarzschild spacetime it is 



Vt 



RW 



1(1 + I) 



6M 



(23) 



and Eq. ([2T|) becomes the standard RW equation. 

One of the main advantages of the GS master scalar is 
that the perturbation of the metric can be algebraically 
reconstructed. But there arc some other variables which 
obey the same RW equation. In particular, one that will 
be very useful for our purposes is the one introduced by 
Regge and Wheeler themselves in their seminal paper . 
Its gauge-invariant generalization to higher orders takes 
the following form, 



(24) 



When using Schwarzschild coordinates for the back- 
ground this definition becomes 



2M -r 
2r2 



f"W-2f"V 



r 



{„} 



0$, (25) 



The second-order source can be written in terms of the 
source of the Einstein equations. 
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^AB \ ^ \ ^ l-e)ci7nmm 
Li fn,m 



(20) 



As in the even parity case, the sources 5* that appear in 
the previous expression are explicitly given in reference 
[l|. Equation ([TO)) is a wave equation for the scalar *"'n. 
It contains all the relevant physical information of the 



where is the standard term that depends on lower 

order perturbations and vanishes when particularizing to 
the RW gauge. There are practical advantages and dis- 
advantages for each of these definitions for the master 
scalar in the odd parity sector, as we will sec in Section 

EH 

At linear order, obeys the same equation as 11. 
But at second order the source term changes. 



(26) 
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IV. RADIATED POWER 



(27) 



plicitly, this will be discussed later) make a gauge trans- 
formation from the RW gauge to an asymptotically flat 
one. Once in that gauge, we will apply Eq. (|28p in order 
to obtain the radiated power. 

It is useful to rewrite the Landau and Lifshitz formula 
()28|) in covariant way with respect to the coordinates on 
S"^, as follows. We first introduce the trace-free projector 
to the sphere (recall that 7q{, is the standard, unit metric 
on the sphere), 



To compute the power radiated to infinity by gravita- 
tional waves we will use the Landau and Lifshitz (LL) 
formula [13, 



dPower 



1 



1 



sm 



dt 



dge 



dt 



'6 dt 



(28) 



This expression is only valid in an asymptotically flat 
(AF) gauge, by which we mean, following |20l |. 



^"^hte-, ^"^ht^, '-"^Ke, ^"X</. 



= 0(r-'), (29) 

O(r-i), (30) 

f">/i^^, f">V = 0{r), (31) 

j-'^-^hab = 0(r°) (32) 



to the desired order n. 

Once the Regge- Wheeler and Zerilli functions are 
known up to a given order, one can reconstruct the per- 
turbed metric components in an AF gauge (we will turn 
to this below). Once those are known, we could in prin- 
ciple use Eq. ([28|) to compute the radiated energy. 

It turns out that the RW gauge is not asymptotically 
flat. Because of that we will need to (explicitly or im- 



Pah''' 
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-slabl 



cd 



(33) 



such that for any rank-two tensor Acd we have that 
Pab^'^cd is trace free. It allows us to introduce the pro- 
jected trace-free metric on the sphere, 

(34) 



— Pab 9cd^ 



so that formula 
way: 

dPower 



can be rewritten in the following 



1 



dt 



jac^bd 



died 
dt 



(35) 



dn 327rr2 
where the star denotes complex conjugation 

Using the perturbative expansion in spherical harmon- 
ics dU, the projected trace- free metric takes the form 



oo 



n=l Lm 



XPab}, (36) 



where the superscript AF stands for any asymptotically 
flat gauge. 

Finally, making use of the normalization shown in ap- 
pendixIXland the fact that the tensor spherical harmonics 
arc trace-free, it is easy to integrate the emitted power 
over the solid angle to get 



J 



Power : 



647rr' 



oo oo 

EE 

3 = 1 k=l 



(^ + 2)! 
j!fc! f-(/-2)! 



E 



dt 



dt 



5 0}/jmAF 

dt 



d f'=>/i™AF 

dt 



(37) 



Therefore, the problem of extracting the radiated 
power with an order of e" reduces to finding the value 
of the time derivative of the harmonic coefficients ***G™ 
and for all fc < n, in an asymptotically fiat gauge. 

Note that in the last formula there is no coupling be- 
tween modes with different harmonic labels. This is so 
because of the integrated character of the total emitted 
power, and the orthogonality between different spherical 
harmonics. This has important consequences when one 
wants to obtain the radiated power up to a given order 
e" consistently, as discussed later. 



FIRST ORDER PERTURBATIONS 



We now turn our attention to the reconstruction of the 
metric components from the first- order master functions, 
and the computation of the radiated energy in terms of 
them. The reason for re-deriving these results here is 
twofold. First, to fix the conventions that we will use in 
the second order case, in which the first order perturba- 
tions will appear as 'source terms'. Second, to sketch in 
a less complicated case the type of calculations presented 
in the next section for second order perturbations. 



6 



After writing down the first order master equations we 
reconstruct the metric in the RW gauge. One can expUc- 
itly transform the result from the latter to an arbitrary 
gauge. In fact, in order to use the Landau and Lifshitz 
formula (|28p we must use an asymptotically flat gauge. 
As we will see, the RW gauge is not asymptotically flat. 
Because of this, following we will first perform an 
explicit asymptotic gauge transformation from the RW 
gauge to an AF one and afterwards apply the LL for- 
mula. 

In addition, following [l^, we use an alternative way of 
computing the radiated energy, which exploits the gauge 
invariant form of the master functions presented in pre- 
vious section, instead of making an explicit asymptotic 
gauge transformation. 

In this section we will remove all harmonic labels, as 
well as the n = 1 labels, since all objects will be of first 
order and correspond to a generic harmonic pair (/, m). 



A. Even parity/polar sector 

At first order, the Zerilli equation is a wave equation 
without sources. 



and the tortoise one r* . The Zerilli function ^ can be ex- 
panded in inverse powers of r, with coefficients depending 
on the retarded time u = t — r* , 



O 



(43) 



The Zerilli equation is then equivalent to a series of rela- 
tions, order by order in r, among those coefficients: 



«'o(u) = 



/(/ + !) 
*2(w) 



F{u), ^,{u) = F{u), 



A 3M(A + 2) . 

— -t (U) ; —I' (U), 

2 ^ ' 2A(A + 1) ^ 



(44) 



(38) 



This master variable contains the physical information of 
the system since it is possible to reconstruct from it all 
components of the perturbation of the metric, in the RW 
gauge. We explicitly display the results of such recon- 
struction: 



2M 



{ 



+ {M'l + Xrf 
2(2M - r)r2 {QM + 2Ar)^ ^" 
+ 4r [A(/2 + Z- 8) r^Af -2AV^ 
+ 4[l8M-'' + 18ArA/2 + 6A2r2A/ + /(/ + l)A''r 

Htt, 



2^31 



{2M 



Htr — 



2 (3A/2 + 3Ar7\/ - Ar^) 



l{l + l)[MP + {P + 1 
2r 



5)rM - Xr- 



1 



2;(/ + l)r2 (3A/ + Ar) 
{ 2r [-12A-/2 -2(1'^ + 1-5) rM + 2\r'^] ^' 



(39) 

(40) 

(41) 
(42) 



24M-' 



12ArM + + 1)(; + 2)r2] ^i] 



Introducing these relations into the linearized Einstein 
equations, one can show that all of them are trivially 
satisfied if the Zerilli equation ([55]) holds. 

Next we explicitly display the divergent nature of these 
quantities (and, as a consequence, of the first order met- 
ric perturbations in the RW gauge). For that purpose 
we temporarily introduce Schwarzschild coordinates (i, r) 



where F{u) = dF{u)/du. The function F{u) can be un- 
derstood as the free data at null infinity. 

In order to see the divergent behaviour of the harmonic 
coefficients in RW gauge at null infinity, we replace the 
expansion ^ in p9ll42p to obtain 



Htt 

Htr 

K 



4F 



4A^^' 



;2(/ + i)2 i-^{i + if 

4'f mM'F+4.XF (I 



4F 



4f 



8M'f +iXF 



O 



(45) 
(46) 
,(47) 
(48) 



where the orders r'^ and r ^ vanish for the harmonic 
coefficient K. 

In order to apply the LL formula we perform an ex- 
plicit asymptotic (that is, near null infinity) transforma- 
tion from the RW gauge to an AF one. We will not show 
the form of the resulting change of coordinates but in- 
stead directly show the asymptotic form of the metric 
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coefficients in the new gauge, 



B. Odd parity /axial sector 



ttAF 



+ O(4l, 



o + o 



AF 



O + O 



$1 



4P(Z + 1)2 



O 



2/2(^ + 1)2 

i ) r 

AF 1 



•• (Z + 2)! 

AMF + - ■ — F 

+ {l-2)\ 



2MF 



/2(/ + l)2r P(/ + l)2r 



1 



4AF 1 „^ 1 
2$i — 



K 



AF 



+ 0(1), 



(49) 
(50) 

(51) 

1 

r 

(52) 



(53) 

(54) 
(55) 



where zeros stand to show that in fact one could ask for 
faster decay rates than the ones defined in (|29ll32p and 
$1 = <i>i(u) is a gauge freedom that it is not fixed by the 
requirement of asymptotic flatness. From the behaviour 
of the harmonic coefficient G in an asymptoticaUy flat 
gauge (|54p and the asymptotic expansion of the ZeriUi 
fimction it is easy to obtain 



G 



AF 



2^- 



l)r V 



(56) 



Alternatively, this last relation can be directly ob- 
tained from the gauge invariant definition of the ZeriUi 
variable (|imi2[) . Since that definition is valid for any 
gauge we can suppose that we are in an AF gauge. Im- 
posing the decay rates (|29ll32p it is straightforward to 



obtain The advantage of this last method is that 

wc do not have to do an explicit asymptotic gauge trans- 
formation. The disadvantage is, however, that we need 
to assume that (|29ll32p is indeed possible. 

Either way, using the relation (j56p and the LL formula 
we obtain the radiated power in terms of the ZeriUi func- 
tion. 



Power : 



l,rn 



2) 



;(/ + !) 



dm 



at 



(57) 



Since this expression holds asymptotically, we can now 
forget that we have used intermediate Schwarzschild-type 
coordinates to derive it, since it will hold for any time 
coordinate which agrees with it at infinity. 



We proceed as in the even parity/polar sector. We first 
reconstruct the metric from the RW scalar 11 satisfying 
the RW equation 



1 



(r^n)! 



2r2 

in the RW gauge: 
ht = 

hr = 



(/-!)(/ + 2) 



n = 0, 



2A 



(2M-r)(4n + rn') 



n. 



2A(2M - r) 



(58) 

(59) 
(60) 



Next we expand the master scalar in inverse powers of 
r near the asymptotic null infinity (r ~* oo,u = const.). 
Since 11 = r^II obeys a standard wave equation, our mas- 
ter scalar will have the following behavior. 



n 



_ no(u) Hi (it) n2(u) 



r.3 



Oil 



fyi'-± r^ij 

We can define a function J(u) such that 
2 



(61) 



Jiu), 



X 



Ui{u) = J{u), 
3M 



2-^^") /(/ + 1) 



Jiu). (62) 



With these expansions at hand, we can obtain the pre- 
cise divergent behaviour of the RW gauge in terms of the 
fimction J(tt), 



ht 

hr 



Xl{l + l) 
r 



J 



l{l 



+ 1 \r 



(63) 



Xl{l + l 



-J - 



As in the even parity/polar sector, we make an explicit 
asymptotic gauge transformation to an AF gauge. And, 
as in that sector, we do not show the details of the re- 
sulting transformation but instead the final asymptotic 
behavior of the metric in the new gauge: 



•AF 



1 



-J+ii^ 

4 (/ + 2)! 
G-2)! 



MJ 



1 



O 



1 



(65) 



O 



AF 



(1 + 2)1 



-J - 



MJ + -{l^ + 1 - . 



(66) 



Xl{l + 1) 1(1 + 1) r 



J+-E0+O f^),(67) 



where Sq = ^o(^i) is a residual gauge freedom. From 
there it is easy to obtain that asymptotically. 



AF 



-n + o(r"). 



(68) 
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Replacing this result in the LL formula for the emitted 
power we obtain 



Power 



IGtt ^ (I- l)(l 



■2) 



dt 



(69) 



One can try to obtain relation ([68|) with a gauge invari- 
ant approach, as we have done in the polar case. It is not 
possible in this sector, however, since the gauge invari- 
ant form of the master variable 11 (fTSl) does not contain 
the harmonic coefficient h. Hence, if using that variable 
one has to go through the explicit gauge transformation. 
But at this point we note that there is another master 
variable <I> whose gauge-invariant form (j25p does contain 
the harmonic coefficient h. Making a transformation to 
outgoing coordinates and assuming that we are in an AF 
gauge we can easily obtain between $ and h at null in- 
finity. 



LAF 



2r4' + C'(r°) 



(70) 



Therefore, the emitted power can also be given in terms 
of this last variable, 



l,m ^ ' 



(71) 



Because we can apply this gauge- invariant approach to 
relate the master variable with the harmonic coefficient 
h at null infinity, at second-order we will use the variable 
^^'<J>. But there is one disadvantage of using <I> instead 
of n. We have shown above the reconstruction of the 
perturbations of the metric in the RW gauge in terms of 
n ([59]) . These relations are algebraic. If we try to do the 
same with the variable <&, we find that the reconstruction 
of the metric is not algebraic, but differential. 



- 2Af ) 
2M 

ht = 1 



hi 



2M hr 
- 2M~ 



(72) 
(73) 



This is why we will use 11 at linear order and <& at second. 



VI. SECOND ORDER PERTURBATIONS 



it obeys, '^^^S'*, diverges at large radii. In order to see 
this, it is sufficient to take its gauge-invariant definition 
PHI12|) . assume an asymptotically flat gauge, and impose 
conditions (|49ll55p and (|65ll67p . In this way, we find that 
the quadratic source ^^'Q* diverges as. 



{2}, 



= Q2r^ + Qi 



oil 

r 



(74) 



where Qq, Qi and Q2 are quadratic functions of 
{F, F, J, J}. The hat and bar on F and J denote the 
generic harmonic labels (l, rh) and (/, m) of two first-order 
modes, respectively. For instance, the dominant term is 
given by 



32 



(75) 



where the iJ-coefficients, defined in Appendix 1X1 are nor- 
malized products of two Clebsch-Gordan factors. In this 
case, the E coefficient restricts the sums to those har- 
monic labels {1,1, 1}, such that l + l + l is an even number. 
That is, for the cases in which Z Hs an odd number 
the term Q2 cancels out. In contrast, the term Qi has a 
non- vanishing contribution in all cases. 

These divergences arc non-physical, and have their ori- 
gin in the freedom to add low-order gauge-invariant terms 
in the definition of the second and higher order Zerilli 
functions, as discussed in Section IlIII Fortunately, they 
can be removed by exploiting that freedom to regular- 
ize the source of the Zerilli equation. In that way 
will obey the same equation but with a different, non di- 
verging source term. The following discussion is rather 
technical, but necessary. The reader not interested in 
the details, though, might skip it and refer to Eq. (jST]) 
[and (inH) for the odd parity sector] , keeping in mind that 
we have made use of the freedom in defining the second 
order master functions in such a way that their associ- 
ated sources are non-divergent both at the horizon and 
at infinity. 

Our aim is to obtain some quadratic terms on the first- 
order Zerilli functions and GS master scalars Qreg = 
Qregil^', ^, II, n}] which rcproducc the asymptotic di- 
vergent behavior of the source near null infinity. That 
is, for r » M with u ~ const, 



In order to solve for the second-order perturbations it 
is in principle enough to solve the Zerilli (fT3)) and RW 
(|26p equations with their corresponding sources [given by 
(fT6|l and ([27|l . respectively]. However, as we will discuss 
in the next two subsections, in practice there are some 
technical obstacles to overcome first. 



A. Sources: even parity sector 

As we have defined it, the second-order Zerilli function 
^^^^^f, and in consequence also the source of the equation 



(76) 



Q,e3[{*,^',n,H}] = Q2[{F,F,J,J}y + 

+Qi[{F, F, J, J}]r + Qo[{F, F, J, J}] + O 



In order to construct the function Qreg we make the fol- 



lowing replacements in Q2, Qi and Qq, 



F 



y^(r 



i)n. (77) 



These rules include all cases but the first and zeroth 
derivatives. There are no or J terms without deriva- 
tives in the divergent terms, but there are some first-order 
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derivatives. Hence, the straightforward definition would 
be 



dr 



(78) 



However, these replacements introduce divergences at the 
horizon r = 2M . In order to see this, we choose ingoing 
Eddington-Finkelstein coordinates, which are smooth at 
the horizon. They are obtained from the Schwarzschild 
coordinates (t, r) by the following transformation 



t^w = t + 2Mln 



2M 



- 1 



(79) 



In these coordinates the two-dimensional background 
metric takes the form 



gAsdx^dx^ 



1 + ]dr^ 

r 



AM 



-dwdr 



(80) 



Therefore, we have the following relation between coor- 
dinate derivatives, 



dr 



9* 

dr 



r + 2M fd-^ 



- 2M \ duj 



(81) 



which makes explicit the divergence of the radial deriva- 
tive in outgoing coordinates at the horizon r — 2M. Tak- 
ing into account Eq. ([78|) this implies that the source 
diverges there as well. 

In order to regularize the source at large radii without 
introducing divergences at the horizon we proceed in the 
following way. First, we make a Taylor expansion in in- 
verse powers of r of the right-hand side of Eq. ((8T|) . Next, 
we define a derivative that approaches {d/dr)u for large 
r, but without being divergent at the horizon. Following 
this method we get 

which is finite at the horizon. Converting this last re- 
lation into Schwarzschild coordinates gives the following 
rules to reconstruct the divergent terms, 



F 



J 



,.2 /,,3n^ 



^ 2M~r ' 
, . r3 - 16Af3 3 . 



2M - r 



-r'^H. 



(83) 
(84) 



The replacements ([77|) and (|83|) must be done system- 
atically. That is, first take Q2r^ and reconstruct the term 
that will reproduce it. 



i,rh l,m 



8r' 



Xl{l + + 1) 



E°%^^SJ. (85) 



When expanding near null infinity, this term will go as 
Q2r^ + Rir + RQ +0{r~-^). In order to remove the diver- 
gent terms of order 0(r), it is not enough to find a term 
that will reproduce Qir, it must reproduce {Qi — Ri)r, 
to compensate the new term we have just introduced. 
Therefore, we take (Qi — Ri)r and make the above re- 
placements ([77)) and ([55]) again. And so on, until we 
achieve the desired quadratic function Qreg[{^, H, H}] 
which asymptotically behaves as in Eq. (j76p . 

In this way, we define the regularized second-order Zer- 
illi function as 



ireg • 



It obeys the following wave equation. 



1^ - Vz '''^re 



m _ {2}greg 



(86) 
(87) 



where the regularized source is given by 



{2]qreg _ {-2} q 



ireg \A 



reg- 



We have implemented this regularization procedure for 
generic (Z > 2) first and second order modes, but the 
results are quite lengthy. Just to illustrate the point, 
though, we explicitly show the final result for the regu- 
larization factor for the particular case {l,rh) = {l,m) = 
(Z,m) = (2,0): 



(89) 



''■'^^ 252 (2Af - r) V TT 

{ 2(2A/ - r) (^{9M + r) + 6 

+ (llOA/^ - 21rM^ + Ur^M + 4r^) ^'^i/ 

- 2(2Af - r) (4r2 - (15M - 6r) f'^*) *''\^| 

/ ^ (i6 {i>ri *^>n + (2r - 3A/) f^>ri <^>ri| . 

: V TT L J 



3r^ /5 
224' 



B. Sources: odd parity sector 

In the axial case there is no such a divergence. Fol- 
lowing the same steps as above, one finds that near null 
infinity the quadratic part of the RW function tends 
to 



{2h 



= Q 



(0) 



O 



(90) 



Therefore, in principle there is no need to regularize the 
second-order RW function. But, as it will be clear in 
the next subsection, we are still interested in removing 
the term of order 0(1). We do so by applying the same 
procedure as in the polar case: namely, we obtain a term 
which reproduces Q^"* at null infinity. 
After that we define the regularized second-order RW 
variable as 



(91) 
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and its evolution equation 

'''^regT\A^^ - ^RW ''''^regT = '^'Sl^' , (92) 

where the regularized source is again given by 

<^^Sl^^ EE '^>5* + Q-^i^l^ - VnwQT''. (93) 

As in the even parity case, we shall not present the 
details of the general procedure but instead simply ex- 
plicitly show the final result for the regularization factor 
for the (Z, m) = {I, rh) = {I, m) = (2, 0) case: 

Q*' = -^\/^{3M + nvi/ + *n}. (94) 

84 V TT 

The regularized sources for the equations of motion 
(|87l) and ([92]) are one of the main results of this arti- 



cle. We have calculated them for the presence of any 
first- and second-order axial or polar modes. We do not 
include their explicit form here because they are quite 
lengthy and do not contribute to the discussion. They 
are available from the authors upon request. As an il- 
lustration, for the interested reader, in Appendix [B] we 
have written down their explicit form for some particular 
values of the harmonic labels. 



C. Radiated power 

Once we have solved for the first and second-order mas- 
ter equations we can obtain the radiated power by using 
Eq. ([37|) . Expanding it explicitly up to order e^, it takes 
the following form 



Power 



(? + 2)! 



647rr2 ^ (I - 2)1 

l,m. 







2 


QiD^niAF 


1 




+ 


+ eRe 




dt 




dt 



dt 



dt 



dt 



lAF - 



dt 



(95) 



where Re means the real part. Again, the problem of 
finding the radiated power reduces to calculating the har- 
monic coefficients G™ and h™ near null infinity, in an 
asymptotically flat gauge. More precisely, we want to re- 
late them with the regularized master scalars constructed 
in the previous two subsections. 

In those subsections we regularized the second-order 
master variables so that the quadratic contributions from 
first-order modes decay as 0{\/r) near null infinity. 
Hence, we can use their gauge-invariant definitions, (jl2[) 
and (fT5|). and assume an AF gauge (|29ll32p up to sec- 



ond order. This leads to the very same relations as at 
first-order; namely, 

{2};^mAF ^ 2rt^'$™re3 +0(r"). (97) 

Replacing these expressions in Eq. ([M]) . the radiated 
power up to order is given in terms of the master 
scalars by 



Power 



^ (/ + 2)! 
647r ^ il- 2)! 

l.ni 

2^.3 afDn 



^2(^ + 1)2 



dt 



A2 



dt 



eRe 



dt 



I reg 



dt 



1_ {2}(j)»" 



A 



dt 



I reg 



(98) 



This last formula, complemented with the evolution 
equations for the regularized master scalars constructed 
above, provides a closed set of formulas that permits to 
obtain the radiated power up to order in the most 
general case in a fully consistent way. 



At this point we want to discuss an aspect of second 
order perturbations of Schwarzschild black holes which 
seems not have been discussed before in the literature. 

Namely, even when we are solving for the metric up 
to second-order perturbations, we can only obtain the 
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complete radiated power up to third-order in e. In order 
to obtain the following order e^, one should also consider 
third-order perturbations. 

Let us further elaborate on this point. Consider the 
simplest possible scenario: a unique first-order mode with 
harmonic labels {l,m) and polarity a. The polarity a 
will take the value 1 for evcn-parity/polar modes and 
— 1 for odd-parity/ axial modes. Because of reality condi- 
tions, if the mode (I, m, a) is present, so is its conjugated 
{l,-m,(T). 

The self-coupling of this mode will generate several 
second-order modes but not in general the one with la- 
bels (l,±TO,cr). In contrast, at third-order the mode 
with indices {I, ±to, a) will indeed be generated. This 
means that the third-order modes will always contribute 
to the emitted power at order e*, coupled to the first- 
order mode with the same harmonic labels. Therefore, 
without considering third-order modes, one can only ob- 



tain the radiated power consistently up to order e^. 

In order for the emitted power (j37p to have a contribu- 
tion of that order (e^) the self-coupling of the first-order 
mode must give a second-order mode with the same la- 
bels (l, m, cr). It is easy to see that, when only considering 
a first-order mode (Z, ±m, cr), this will happen if and only 
if m = and if, for cr = 1 (cr = —1), Z is an even (odd) 
number. 

In order to make the above discussion more explicit and 
analyze which problems can be addressed consistently, 
let us consider the particular case of a first-order even- 
parity/polar mode with harmonic labels I = m = 2 and 
I = m = —2. These modes will generate the second-order 
{Z = 4,m = ±4,0}, {Z = 2,m = 0} and {I = 0,m = 0} 
even-parity/polar modes as well as the {Z = 3,to = 0} 
odd-parity/ axial mode. Particularizing the power for- 
mula (|37[) in terms of the master scalars to this case, we 
obtain the following contributions from the modes. 



J 



Power = |i9t*'^^'2l 

127r I ^1 



9£' 



640 



r 



15£^ 



96^ 



(99) 



where the second-order master scalars are the regularized 
ones. Here it can be clearly seen that the order e'^ is not 
present. The problem with this last formula is that it is 
not complete since the third-order {I = 2, m = ±2} polar 
mode would contribute to the power at order e^. 

On the other hand, let us consider the first-order 
mode / = 2 with all its possible harmonic labels m = 



0, ±1,±2. By coupling, they will generate the second- 
order polar modes / = 0, / = 2 and / ~ 4 with all 
their possible m. That is, we will have the second- 
order {Z = 0,TO = 0}, {I = 2,m = 0,±1,±2} and 
{Z 4, m = 0, ±1, ±2, ±3, ±4} polar modes. This partic- 
ular case will provide a non-vanishing e^-order term to 
the power. 



J 



Power = ^|2|at<^>*2|2^2|9t<^>^i|'-f|at*^>*°|'} 

+ ^Re [2dti'''^l)dti'''^lr + 2dt{'''^l)dt{'''^lr 



9t(*^'*^)at(f^%0)* 



(100) 



r 



where, again, the second-order Zerilli function must be 
understood as regularized. In this last case the formula 
is exact up to the displayed order, to which the generated 
second-order axial modes and third-order polar modes do 
not contribute. 



VII. FINAL REMARKS 

In this paper we have introduced a complete gauge in- 
variant formalism to study arbitrary perturbations of a 
Schwarzschild black hole up to second order. In particu- 



lar, we regularized the resulting equations, making them 
suitable for a numerical implementation. 

This formalism enables a variety of applications and 
studies. These range from the non-linear stability of the 
black hole horizon, to non-linear features in gravitational 
waves and mode-mode coupling. We will report those 
studies elsewhere. 

All calculations of this paper have been done with -and 
have been largely possible at all due to- new, efficient 
symbolic manipulations tools. The resulting expressions 
in most cases are very long and their explicit expressions 
are not particularly enlightening. For that reason we 
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have refrained from explicitly presenting most of them. 
They are however, available upon request. 
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(STF) tensors 



STF 



X 



I ai...as 



{Yi™:ai...aJ 
b ym 



(Al) 
(A2) 



together with the metric 7o& and the antisymmetric ten- 
sor €ab dl- For the particular case s = 0, those objects 
must be read as ee and XJ^ = 0. They are 
normalized in the following way. 



dnxrabY'^i'"' [xjv'cd 



I {I + 2) 



2 {I -2) 
I {I + 2) 



2 {I -2) 
0. 



(A5) 



Going beyond linear perturbation theory, the nonlinear 
coupling between two first-order modes results in prod- 
ucts between two tensor spherical harmonics (Z, m, s) and 
[V , m', s'). Those products can be decomposed into a lin- 
ear combination of harmonics {I" ^m-\-m' ,s + s') with an 
explicit formula involving coefficients 



^s'l'TTl'l" 



k{i',\s'm,\s\) 



C, 



m' mm -{-m ^s' s s' +s 



I' I I" 



k{l",\s + s'\) 

arc the usual 
and fc is a normalization factor defined by. 



1 1'' 



(A6) 



where C™ arc the usual Clebsch-Gordan coefficients 



k{l,s) 



' {21 + + sy. 
2^+2 7r(/- s)!' 



(A7) 



These ^^-coefficients encode the geometric selection rules 
that determine which pairs of modes do actually couple. 
See [H for fuU details. 



APPENDIX B: REGULARIZED SOURCES, SOME 
EXPLICIT EXAMPLES 



APPENDIX A: TENSOR SPHERICAL 
HARMONICS 

Tensor fields of any rank s on the sphere will be decom- 
posed using a basis of tensor spherical harmonics. Such 
basis can be constructed from the symmetric trace-free 



In this appendix we show two particular examples for 
the regularized sources of the RW (|92p and ZeriUi equa- 
tions ((871) . 

Let us first assume that we have the first-order {I = 
2,777 = 1} cven-parity/polar and {I = 8,m ^ —4} odd- 
parity/axial modes. The regularized source generated 
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by them for the equation of motion of the, for example, given by 
second-order {/ = 7, m = 3} even-parity/polar mode is 



+ 6on,„ (e;/^ + bbu^ + 7U - isf r^-^^tr 

- 20n triU + 9f {18U^ + 51C/3 - - 26U + 20) r^^^r 

- 5n,t(3C/ + 2)2 (212C/^ + 3364C/^ + 11603C/2 - 13149C/ -f 3240) r^-^^rr 
+ 20n,„.(C/ + 9)2 (126C/** + 141C/^ - 82U^ - 50U + 20) r^'^^t 

+ 5n,r(3[/ + 2)2 (308C/^ + 4972C/^ 17255C/2 - 22221f7 + 6156) r^^-^^^ 

- 180n tr {U + 9)2 (6C/^ + 9U^ + 2C/2 + 4[/ _ 4) 

+ n t {360U^ + 2568U^ + 57529?/'* - 14036;73 + 375894[/2 + 88254C/ - 123444) r^',^ 

(14220C/'^ + 253302C/^ -f 1234181C/'* + 854111C/^ - 966354[/2 - 375354C/ + 220644) r^^t 

+ 15n(3C/ + 2)2 (92f7^ -f 1492^73 + 8507C/2 + 14154C/ - 9396) r*,*^ 

- n,t (6840C/^ + lUUrn^ + 422069c/* + 1424C/^ - 213006C/2 + 271944C/ - 48924) ^ 

+ n (6210C/^ + 116313C/^ -f 1283789C/* + 6929894C/^ + 6649074C/2 - 316926C/ - 1359504) *,t|, (Bl) 



where the symbol U stands for the dimensionless mass 
U = M/r. 

As a second example, we show the regularized source 
for the RW equation for the particular case in which 



the first-order even-parity/polar modes (/ = 3,m = 0) 
and (/ = 4,m = —1) generate a second-order odd- 
parity/axial mode with labels (/ = 4,to = —1): 



S. 



reg 



3i 



I - lOr (3[/2 + UU + 15)^ ^,,.„*,^^(2C/ - 1)^ 



8800V7^([/ + 3)4(2[/ - 1)2(3[/ + 5)^ 
+ 26r (3C/2 + 14C/ + 15)"^ *,^^*,„.r(2[/ - 1)^ 

~ ^^^^ (3060C/^ 4- 49401C/^ -f 332356[/6 -f 1197973C/^ + 2636572C/'' 
-h 3760905C/3 + 2764530C/2 - 467775C/ - 1518750) ^,„§(2C/ - 1)^ 
+ ^^^2^^ (1620C/** + 28161C/^ + 173844C/6 + 197637[/^ - 1511900[/* 
- 5534775C/3 - 7023550C/2 - 3510375C/ - 573750) **,^^(2C/ - 1)^ 

+ lOr (3L/2 + 14C/ + 15)"* *,t„*,tr(2?7 - 1)^ - 26r (3C/2 + 14J7 + 15) if ,tr^ ,trr{2U - 1) 
16(1 - 2C/)2 



;;i701C/" + 35262t/^" -f 320166C/^ + 1720086C/* + 6285736C/^ 

16821825C/*' -t- 34748135C7^ + 56990175C/* 4- 68601150C/^ + 42931125t/2 
5703750t/- 15946875)** 

^ (6[/2 + 7L/ - 5)^ (1530C/^ + 34221t/^ + 303099?/^ -h 1485635/7*^ + 4592169C/'^ 
9179205L/-* + 10353033[/3 + 3316365[/2 - 2994975C/ - 1478250) ^.^^ 

10 4 * _ 

— ([/ -I- 3)2 (6[/2 7C/ - 5) 9C/2 27C/ + 90) *,r,.r* 
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+ ^ + 5C/ - 3)^ (24138C/^ + 399357C/*^ + 2535795?/^ + 88662631/^ 

+ 20189321C/^ + 31979265J7'' + 34936825C/3 + 20024625C/2 - 4674375C/ - 9618750) 

+ ^ (6[/3 + 25C/2 + 16U - 15)^ (90C/^ + 468f7^ + 1763C/^ + 5632C/4 + 221Q4U^ 

- 6480[/2 _ 35025[/ + 15750) ^f ,,^^ 

+ ^(C/ + 3)2 (6[/2 + 7?/ - 5)^ (675C/^ + 5741C/'^ + 15946C/^ + 24570C/2 + 10590C/ - 13950) ,rr^ ,r 

- 5(C/ + if{llU - 9) (6C/2 + 7C/ - 5)** ,rrr^ ,r 

- ^(3C/ + 5)2 (2[/2 + 5C/ - 3)^ (1377;/^ + 11403C/^ + 25994C/3 + 19250C/2 _ ^iqjj _ ^^0350) ^-^^^-^^^ 

- 8(1 - 2Uf (3t/2 + UU + 15)^ (30C/^ + 91C/2 + 99C/ + 15) ii ,rr^ 

78 4 - 

+ y(3;7 + 5)2 (2[/2 + 5C/-3) (3[/^ + 15;72 + 25t/ + 50) «'^',rrr 

+ 13(3C/ + ^f{nU - 15) (2[/2 + 5[/ - 3)*^ ^^r'^^.rrr 

- -| (3[/2 + 14C/ + 15)^ (198[/^ + 1164?/^ + 2627C/'^ + 17137?/^ + 45972C/3 
+ 3140C/2-38100C/+ 11475) 'I',**,* 

- ^(C/ + 3)2(3C/ + 5)3 (474?/^ + 4575?/^ + 13106[/* + 20156C/3 + lOUt/^ _ 23685?/ + 8IOO) *,tr*,t 
+ 5(1 - 2UY{U + 3)3(3;/ + 5)^ (4^2 ^ 23^ - 9) ^.t^^*,* 

+ ^([7 + 3)3(3C/ + 5)2 (2142C/^ + 130G5C/^ + 21826C/* + 13968C/3 - 437GC/2 - 21115C/ + 8IOO) ^.t'^.tr 

+ 24 (^2 _ 3[/ _ 5) (6[/3 ^ 25C/2 _^ -^gf; _ -^5-)3 

- 13(l-2C/)2(t/ + 3)''(3t/ + 5)3 (12C/2 + 37C/- 15)*,t^,t^,r|. (B2) 

I 



